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ABSTRACT 

We show that instantons violate a supersymmetric identity in a clas- 
sically supersymmetric Higgs model with no massless fermions. This 
anomalous breaking arises because the correct perturbative expansion in 
the instanton sector is not supersymmetric. The attempt to construct a 
manifestly supersymmetric expansion generates infra-red divergences. 
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1. Introduction 



Asymptotically free supersymmetric gauge theories may play an important role 
at very short distance scales. Present day understanding of the physical properties 
of these theories [1-6] relies in a crucial way on the conjecture that non-perturbative 
effects do not break supersymmetry (SUSY) explicitly. The plausibility of this conjec- 
ture has been discussed in the literature. On the other hand, there are also arguments 
against it 0, [§] • 

In this paper we show that the one instanton contribution to a certain condensate 
violates the relevant supersymmetric identity. We consider a SUSY-Higgs model, 
which allows for a fully controlled calculation of non-perturbative effects. The model 
contains no massless fermions. Hence, the option of spontaneous "dynamical" SUSY 
breaking is ruled out. 

Anomalous SUSY breaking is closely related to the Higgs mechanism, which pro- 
vides a physical infra-red cutoff through the generation of masses. Physical observ- 
ables receive non-supersymmetric contributions which arise from the non-constancy 
of the Higgs field in the instanton sector. 

A formally supersymmetric expansion can be constructed in the instanton sector 
if the effects of the classical Higgs field, in particular the generation of masses, are 
incorporated perturbatively. Only the background gauge field is taken into account at 
tree level. In this formal expansion one should find supersymmetric results provided 
the results are finite. 

The price paid for the perturbative treatment of mass terms is the appearance 
of infra-red (IR) divergences. In the literature on supersymmetric QCD there are 
several leading order instanton calculations with a non- vanishing scalar VEV || ||. 
In all of them, IR divergences will arise at the next-to-leading order. In the case of the 
condensate calculated below, there are IR divergences already at the leading order. 
Moreover, the IR divergences become more and more severe as the order increases. 

The anomalous breaking of SUSY is triggered by these IR divergences. While 
the supersymmetric calculus is valid for some (but not all) leading order calculations, 
the supersymmetric Feynman rules are in general ill-defined. As a result, an IR 
cutoff must be introduced explicitly. But the IR cutoff necessarily breaks SUSY if 
we are interested in the one instanton sector]]. This sets the stage for an anomaly 
of a new kind. As we discuss in detail below, the IR divergences are eliminated by 

1 We comment that periodic boundary conditions or the compactification of R 4 to S 4 are both 
incompatible with the asymptotic behaviour of the Higgs field, which reflects the topology of the 
instanton sector. 
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a resummation procedure that reconstructs massive propagators from massless ones. 
This resummation gives rise in general to non-supersymmetric results, because the 
spectrum of massive modes in the instanton sector is not supersymmetric. 

Let us describe the massless supersymmetric expansion in more detail. As the 
background gauge field one takes the pure instanton expression. (The source current 
in the gauge field's equation of motion is neglected to leading order). For simplicity 
we will henceforth assume that all masses arise through the Higgs mechanism. Mass 
insertions are treated perturbatively. Thus, zero modes and propagators are obtained 
by solving the massless field equations in the pure instanton background. The classical 
Higgs field itself is a solution of the equation D 2 (p = 0, which satisfies the boundary 
condition (j){x) — > v at infinity. In this framework, supersymmetric relations exist at 
two levels. As shown in ref. [||] the classical fields and the fermionic zero modes are 
related via SUSY transformations. Also, supersymmetric relations exist between the 
propagators, because the continuous spectrum in the pure instanton background is 
universal 0. 

The supersymmetric expansion reflects an approximation which is valid in the 
instanton's core. In a Higgs model, the leading order effective action in the instanton 
sector is || 

S E (p)=8rc 2 /g 2 (p)+4rc 2 v 2 p 2 . (1) 

Here p is the instanton's size. The numerical coefficient of the second term, which 
represents the contribution of the Higgs field(s), is valid for an SU(2) gauge group [[|. 

Eq. ([[]) implies that the saddle point of the p- integration occurs at p ~ v" 1 . In 
units of the Higgs VEV, the gauge field's strength in the core is 0(l/g). In contrast, 
the magnitude of the Higgs field is 0(1), namely it is comparable to a typical quantum 
fluctuation on the same scale. This provides a formal justification for keeping only 
the gauge field at tree level. Propagators, insertions of the classical Higgs field and 
the full set of zero modes of the pure instanton background are all treated as elements 
for the construction of Feynman graphs. 

The IR divergences arise from multiple mass insertions on a given line in a Feyn- 
man graph. Their physical origin is easily understood. Consider the linearized field 
equations taking into account all classical fields in the usual way. These differential 
equations determine the exact form of propagators and zero modes. In the core, 
the equations are dominated by the covariant derivative terms. Consequently, the 
supersymmetric expressions provide good first approximations for the exact ones. 

In the transition region v~ l <C \x\ <C mT 1 the equations are dominated by the 
free kinetic term. (Here m stands for a generic Higgs- induced mass). However, for 
\x\ ~ m _1 and at larger distances, the mass terms are comparable to the kinetic 
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term. At any distance scale which is large compared to the instanton's size, the 
massless field equations give rise to a power law behaviour. But since the relevant 
fields are massive the correct asymptotic behaviour is a falling exponential. Within 
the supersymmetric Feynman rules, one must sum over an arbitrary number of mass 
insertions on every line in order to reproduce the falling exponential. The sum of a 
finite number of mass insertions will in general be IR divergent. 

The IR divergences of the massless supersymmetric expansion are therefore spu- 
rious, in the sense that they arise from using expressions that do not have the correct 
asymptotic behaviour. This statement applies to all elements of the perturbative ex- 
pansion, namely, to propagators, to zero modes and to the deviations of the classical 
fields themselves from their limiting values. 

The massive Feynman rules in the instanton sector are obtained by applying the 
standard variational principle. To make contact with the massless supersymmetric 
expansion, it is convenient to separate out the integration over the instanton's size. 



This is done by first fixing the size by introducing a constraint [K| . For any value of 
the constraint parameter one sets up a perturbative expansion in a standard manner. 
In the end, one integrates over the instanton's size which is now represented by the 
constraint parameter p. A judicious choice of the constraint will modify mainly the 
classical gauge field's equation. Physically, the source current arising from the matter 
part of the lagrangian tends to shrink the instanton's size, and this effect is balanced 
by the contribution of the constraint. 

The classical fields obtained via this procedure tend exponentially to their vacuum 
values at large distances. The exponential asymptotic behaviour also characterizes 
the massive propagators and the surviving exact zero modes. As in scattering theory, 
this allows the removal of the finite volume cutoff, and one can work directly in the 
infinite volume limit. 

The massive Feynman rules feature a supersymmetric set of vertices, because the 
classical lagrangian is supersymmetric. But in general there are no supersymmetric 
relations between the classical fields and the zero modes' wave functions^. Moreover, 
the massive fluctuations spectrum is no longer supersymmetric (except approximately 

2 Models with Nf — 1 are a special case. An example is the simpler Higgs model obtained by 
dropping the "lepton" families from the model of Sect. 2. In this case, supersymmetric relations 
between the classical fields and the zero modes survive to some extent also within the massive 
Feynman rules. This is a reminiscent of the situation in less than four dimensions, where the 
topologically non-trivial objects are exact solutions of the classical field equations. In both cases the 
SUSY violation arises from the absence of local supersymmetric relations betweens the spectra of 
bosonic and fermionic fluctuations. As a result, however, there are no SUSY violations at the level 
of tree diagrams. See ref. JjJ for more details. 
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at very short distances). The linearized bosonic field equation and the square of the 
Dirac operator define two Schrodinger-like operators. Due to the non-constancy of 
the Higgs field in the instanton's core, these Schrodinger operators involve different 
potentials. Consequently, massive bosonic and fermionic eigenstates are not related 
by the action of a local differential operator [0. Therefore, the breaking of SUSY is 
ultimately attributed to the non-invariance of the path integral measure. 

We summarize the relation between the massive Feynman rules and the massless 
ones by considering the example of an exact zero mode, i.e. a zero mode that survives 
the introduction of the Higgs field. Other elements of the Feynman rules (propagators 
etc.) exhibit analogous behaviour. The l.h.s. of any field equation can be written as 
Hq + V, where V contains the dependence on the Higgs field and Ho contains all the 
rest. Let \l/o be a zero mode of H and \I/ the corresponding zero mode of H + V. 
Formally, V is small compared to H , and so we may think of reexpanding \1/ using 
the Born series ^ = — GqV^/o + • • •■ But in general this is not possible because 
V contains a mass term that changes the asymptotic behaviour. Individual terms in 
the Born series have divergent norms, and their insertion in a Feynman graph gives 
rise to IR divergences. 

In practice, one has to distinguish between two cases. Suppose that the massless 
Feynman rules give rise to a finite answer for some observable (or, more generally, for 
some sub-diagram) at the leading order. As far as a given zero mode is concerned, 
this will typically happen when the leading order result involves only the "zeroth 
approximation" ty which is obviously normalizable. Notice that the difference — \& 
is normalizable too. Replacing \l/o by \l/ — \l/o will measure the difference between the 
correct result and the (still finite) prediction of the massless expansion. One can show 
that, up to logarithmic corrections, the contribution of the difference \1/ — \l/ will be 
damped by two powers of the coupling constant (s). A similar pattern is found for 
other elements of the Feynman graphs. In conclusion, if the massless supersymmetric 
rules give rise to a finite answer for a (sub-)graph at the leading order, then the answer 
is valid at that order, and it will be reproduced by the massive Feynman rules. 

Whenever the massless supersymmetric rules give rise to IR-divergent expressions, 
the massive Feynman rules must be invoked. The mass terms then provide a physical 
cutoff for the formally IR-divergent integrals. Below we encounter a sub-graph which 
is equal, within the massless supersymmetric rules, to the norm squared of the first 
Born correction \&i = —GqV^o to a certain zero mode. However, the norm of is 
IR-divergent. 

Physically, represents a new field component that the original zero mode 
develops due to the mass terms. Using the massive Feynman rules, what has to be 
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calculated is the norm of the new component. (The overall normalization of the zero 
mode is still determined by the original component). The new component's norm has 
a leading logarithmic piece which is easily computed, and which contributes to the 
final result eq. (|]). 

The IR divergences signal the onset of a non- analytic dependence on Higgs masses, 
and, hence, on coupling constants. The condensate calculated below features a non- 
analytic dependence on two Yukawa couplings. This result is of course beyond the 
scope of the massless supersymmetric expansion, which contains only positive pow- 
ers of the Yukawa couplings. More generally, we expect the appearance of analo- 
gous non-analytic dependence on the gauge coupling too. Technically, disentangling 
non-analytic short distance effects (the RG flow) from long distances ones is more 
complicated in this case. 

This paper is organized as follows. In Sect. 2 we define the SUSY-Higgs model. 
In Sect. 3 we present our results. Sect. 4 contains our conclusions. Some technical 
detail are relegated to an Appendix. 

2. The model 

The charged fields of our SU(2)-Higgs model are the same as in supersymmetric 
QCD with N c = Nf = 2. Two charged doublets play the role of Higgs superfields. 
The other pair of charged superfields together with a pair of neutral ones make two 
supersymmetric "lepton" families. The Higgs superpotential is 

Wt = h<$> Qe^eAB ®iA®jB - v 2 ^j . (2) 

Here $ — (^o^o) is a neutral chiral superfield. = (<fiiA, ipiA) contain the two 
Higgs doublets and their fermionic partners. The indices A, i = 1,2 correspond to 
SU(2) colour and flavour groups respectively. (The flavour SU(2) plays the role of 
SU(2)# in the analogy to the supersymmetric Standard Model). We use the following 
representations T AB = — \o~ BA and = for the colour and flavour generators 
respectively. 

The two extra charged superfield are denoted t]a±, and the two extra neutral 
ones are These letters will also be used to denote the fermionic components. The 
scalar components are denoted fjjt± and £i±. The ± index labels the two "lepton" 
families. Notice that these families form a doublet under a horizontal SU(2). The 
full superpotential is 

W = W l + W 2 , (3) 
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where 



W 2 = y e^AB $jB (&+ Va- - &- Va+) 

+ m eij • (4) 

The mass parameter m will be treated as a small perturbation^ and we will work to 
first order in itlq. 

The classical potential has a unique supersymmetric minimum (up to colour and 
flavour transformations). The only non- vanishing VEV is ((pi a) = vSiA- This min- 
imum breaks the gauge symmetry completely, but it leaves unbroken the diagonal 
SU(2) generated by T a + F a . Under this vector SU(2), the two Higgs superfields 
decompose into a sing let = 6iA$u/V2 and a triplet $ a = a? A $ iA /V2. 

All fields acquire masses through the Higgs mechanism. In the triplet sector 
(which includes the gauge and the $ a supermultiplets) the mass is fi — gv. The mass 
of the singlet fields $' and $ is m = \r2hv. For m = 0, the mass of the lepton 
families is m\ = yv. 

The model has two approximate i?-symmetries. U(1)r is a non-anomalous sym- 
metry, which becomes exact in the limit m = 0. The other one, denoted U(l)x, is 
a classical symmetry of the full lagrangian, but it is broken explicitly by instantons. 
The fermion charges under the i?-symmetries are given in Table 1 in units of the 
gaugino's charge. As usual, the charges of the corresponding scalars are related by 
Q n(scalar) = Q #(/ 'ermion) + 1. 





Aa Va± &± 


u(i)« 


-11-1 1 


U(l)x 


-110 



Table 1: Fermion charges under the i?-symmetries 



3 The perturbative treatment of mo does not generate any IR divergences because all the fields 
are already massive for mo = 0. Our results can also be regarded as a calculation of the integrated 
form of the r.h.s. of eq. (^) in the model with mo = 0. 
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3. Results 



3.1 The SUSY identity 

Let 0{x) denote the lowest component of a gauge invariant chiral superfield and 
let Q be a SUSY generator. Using [Q, 0(x)] = one arrives at the following on-shell 
SUSY identity 



= {0{z) {Q, C&-.(x)}] 

= (O(z) ((1/2) ^ + ^(x)+ye jk e AB a + <p kA fj B+ (x))) . (5) 

This identity is related to the analytic properties of the condensate (O) = (O(z)) as 
follows. Let us assume that mo in eq. (§) is a complex parameter. On the second row 
of eq. (HP we contract with and integrate over x. A similar identity is obtained by 
the interchange of £ + <-> and r] + <-> —r\_. Taking into account the two identities 
we arrive at the variation of (O) with respect to m* Q . Therefore, eq. (||) requires that 
(O) be a function of mo only, but not of mj ||]. 

3.2 Calculation of (0 O ) 

Our main result is that the SUSY identity eq. @ is violated for O = 0o- An 
explicit calculation gives rise to 

^=<v^hh logy+ ---- (6) 

The dots stand for subleading termsQ. A is the one-loop RG invariant scale of the 
theory. We expect similar violations for other condensates, such as for example (AA). 
However, for reasons that we explain below the calculation of the mQ-dependence of 
(AA) is more complicated. 

In the calculation of (0o) we adopt the following strategy. Diagrams will first be 
drawn using the massless supersymmetric rules, and we will try to apply these rules to 
evaluate every sub-graph. In view of the relation between the massive Feynman rules 
and the massless ones, as discussed in the introduction, the result for a given sub- 
graph is valid as long as it is finite. When the massless rules produce IR divergent 
expressions, the massive Feynman rules will be invoked to calculate the relevant 
integrals. 

The diagrams that contribute to the leading order result eq. (|J) are depicted in 

Figs. 1 and 2. The two diagrams correspond to the two terms on the last row of 

4 The anti-instanton sector makes a contribution to (<Pq) that satisfies the on-shell relation ((/>g) = 
(<f>o) . We comment that (</>o) vanishes to all orders in perturbation theory because 4>o carries a non- 
zero V(l)x charge. 



eq. (H). The vertex marked with a thick cross is linear in mj. The dotted lines that 
end with a cross represent insertions of the classical Higgs field. They differ from 
ordinary mass insertions by the presence of short range potential terms that arise 
from the non-constancy of the Higgs field in the instanton's core. 

Each thick line that emanates from the instanton (the shaded circle) represents 
one of the eight zero modes that the model had had in the absence of the Higgs 
field. These include the four gaugino zero modes and one zero mode for each charged 
doublet. With the classical Higgs field turned on, only the X s s pair and the T)± pair 
remain exact zero modes (with modified wave functions that contain new field com- 
ponents). The Higgsino and X sc zero modes mix through the Higgs field. However, 
they are still approximate zero modes because the Higgs field is a small perturbation 
over their support. More precisely, the X sc and Higgsino zero modes become either 
resonances or bound states with non-zero eigenvalues. (The former possibility is more 
likely). This is the physical justification for treating their mixing using the massless 
Feynman rules, with the usual reservation that this does not give rise to IR diver- 
gences. Further details on both the exact and the approximate zero modes are given 
in the Appendix. 

Before we continue with the calculation of Figs. 1 and 2 let us discuss the role 
of other diagrams. If we only count powers of coupling constants that come from 
the vertices and ignore the IR divergences momentarily, we should consider diagrams 
which are 0(g A y 2 h). Now, within the massless Feynman rules, all these diagrams 
have in common a quadratic IR divergence due to the external 0o leg- If we try to 
use a (free) massless propagator instead of the correct 4>o propagator, we encounter at 
some point an integral of the form / d A x/x 2 . As we discuss below, within the massive 
Feynman rules the quadratic divergence boils down to a factor of 1/m 2 . We then use 
m 2 = 2h 2 v 2 to separate out a dimensionful factor of v 2 . The remaining 1/h 2 , together 
with the explicit factor of h coming from the vertex, gives rise to the 1/h factor in 
eq. ©. 

Apart from the common quadratic IR divergence, there are several diagrams that 
contain an extra logarithmic IR divergence of the form / d A x/ (x 2 + p 2 ) 2 . The physical 
IR cutoff of such an integral is an inverse mass scale. Thus, the integral behaves like 
the logarithm of the relevant Higgs induced mass times p. In view of the relation 
p ~ v~ 1 , this is equivalent to the logarithm of a coupling constant. 

An enhancement factor of logy arises if and only if the logarithmic divergence 
comes from the "lepton" sector, the relevant diagrams being Figs. 1 and 2. As our 
explicit calculation shows, the leading logs do not cancel between the two diagrams, 
which sum up to the final result eq. (@). 
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Notice that we have indicated in the figures which insertion of the Higgs field 
corresponds to (piA and which to its complex conjugate. Other diagrams that are 
formally of the same order can be obtained by replacing one insertion of the Higgs 
field and one insertion of its complex conjugate by a (00*) propagator. However, 
the resulting diagrams no longer have the extra logarithmic divergence, because the 
propagators provide extra powers of 1/x 2 (some examples are discussed in part III of 
the Appendix). The leading logs arise only from diagrams with a maximal number 
of insertions of the Higgs field. 

Within the massless Feynman rules, there are also individual diagrams containing 
logarithmic IR divergences related to the X sc and Higgsino zero modes. But now the 
offensive l/(x 2 + p 2 ) 2 behaviour cancels out between different diagrams before the p 
integration. This implies that the corresponding diagrams of the massive Feynman 
rules contain no logarithmic factors^ The physical reason for this behaviour is the 
following. A slow 1/x 2 decrease of a zero mode's wave function at intermediate 
distances can arise only from new field components not present in the original zero 
mode. The new subleading field components are well defined in the case of an exact 
zero mode. But subleading field components of approximate zero modes cannot be 
disentangled from the continuous spectrum. For further details see part III of the 
Appendix. 

We now return to the derivation of eq. (^). The diagrams are evaluated with the 
measure (see e.g. ref. @]) 

2VAV 8 / dpV e" 4 ^ V ... (7) 
For Fig. 1 the integrand of the p-integral is 

(ml(yvpf\ogy) . (8) 



'g 2 v 2 \ ( g 2 p 2 \ ( m , N2 . 



V 2 / V 8h , 

Let us explain this expression. Consider first the mixing of the X sc and Higgsino zero 

modes || through the Yukawa-gauge coupling ig\/2 \ a (fi A T A l B ipB ■ To leading order 

there are no IR divergences, and one can calculate the mixing using formal first order 

perturbation theory. Substituting the pure instanton expressions for the zero modes 

eqs. ( |A.6| ) and ( |A.7| ), and for the Higgs field eqs. ( |A.2| ) and ( |A.4| ), we arrive at the 

first factor in eq. (|8|). 

The second factor corresponds to the 0o tadpole, namely to the 0o line that 

emanates from the pair of X ss zero modes (see Fig. 3). It is convenient to trade the 

5 Had the logarithmic divergences coming from the gaugino-Higgsino sector not cancelled out, 
they would have represented an independent contribution to (tfio), because they could not give rise 
to logy, but only to a linear combination of logg and log/i. 
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integration over the instanton's collective coordinates with an integration over the 
external point z, keeping the instanton at the origin. The tadpole is 



h 



-tkitijtAB J d xd zG(z,x)ip iA (x)ip jB (x) . (9) 

Here G(z, x) is the 0o propagator and ^^(x) is a new field component of the k-th X ss 
zero mode (k = 1,2). We immediately see that the attempt to treat 4>o as a massless 
field gives rise to a quadratic IR divergence. Within the massive Feynman rules, the 
0o propagator is defined by 

(-D z + m 2 + U)G(z,x) = 5 A (z-x). (10) 

The short range potential U is a function of the radial coordinate r = \z\ only. 
Explicitly U(r) = m 2 ((p 2 (r) — 1). (See eq. (|A.2 ) for the definition of <p(r)). The 



short range potential can be neglected to leading order. We therefore substitute a 
free massive propagator for G(z,x) and find 



d*zG(z,x) = (11) 
m z 



Notice that the result is independent of x. 

The rest of eq. (|9]) does not contain any divergent factors. Instead of using the 
exact wave function of the ijj(x) component of the X ss zero mode, we can apply first 
order massless perturbation theory. The general formula for the first Born correction 



is \&i = — Go^^o- In the present case Go stands for the (4>4>/ Q propagator in the 
pure instanton background, V stands for the Yukawa-gauge coupling as a function of 
the Higgs field, \l/o stands for the original X ss zero mode eq. (|A.5|) and \?i stands for 
the new i/j(x) field component. In practice, instead of finding ip(x) by integration, it 
is easier to solve the corresponding differential equation H ^i = — V^o, where H is 
the massless covariant Dirac operator. The result is 

#*.(.*) = ¥7T*4r- (12) 

Here a is the spinor index. Carrying out the ^-integration and using m? = 2h 2 v 2 we 
arrive at the second factor in eq. @. 

The last factor in eq. @ corresponds to the upper part of Fig. 1. It takes the 
form 

m* eij J d 4 x £ i+ (x)£j-(x) , (13) 

where £±(x) is a new field component of the 7] T zero mode. As before, we try to use 
the first Born correction. The result is 

= tItHt^. ( 14 ) 

2tt x z + p z 
11 



In this case the spacetime integration is logarithmically IR divergent. The true form 
of £±(x) is determined by eq. (A. 12). It agrees with eq. (|TJj) for |x| <C m^ 1 , whereas 
for \x\ nil 1 , £±(x) becomes a falling exponential. (That the transition occurs at 
the distance scale mj -1 is a general feature of massive linear equations. See part II of 
the Appendix for more details). 

Let us write d A x = r 3 dr dfl, and split the radial integration in eq. ( p^3| ) into three 
regions 

+ J (15) 

The integral is dominated by the intermediate region, where one has £ + £_(x) ~ 1/a; 4 , 
leading to 

H d 4 x 



— j- = —2tt log(mxp) . (16) 



x 4 

Notice that log(mip) = logy+\og(vp), and that \og(vp) is 0(1). Neglecting 0(1) con- 
tributions, we complete the calculation by adding the appropriate prefactors required 
by eq. (fH]). The result is the last factor in eq. (||). 

The calculation of Fig. 2 is similar. Using massless perturbation theory we find 
the "induced scalar field" 

f,Ux) = - % eklX ^ :i . (17) 

This induced field is obtained by solving the massless field equation with a source 
given by the product of the fc-th \ sc zero mode and an 77 zero mode. The second 
induced field, whose source is the product of the n-th Higgsino zero mode (n = 1, 2) 
and the other rj zero mode, is 

Now we have a logarithmic IR divergence from the integral 

m*ye AB f d A x &l(x)(/) iA (x)rj% ± (x) . (19) 



As before, the leading log is obtained by substituting eqs. (|T7D and (|T8|) in eq. ( |19D 
and imposing a physical cutoff at r ~ mj" 1 . Summing the leading logs from the two 
diagrams we find 

(0o) = mSAVi^logy / dp" p^^^p 2 - 1) . (20) 

In this equation, the first and second terms in parenthesis correspond to Figs. 1 and 2 
respectively. Carrying out the p-integration we finally arrive at eq. (pi). Further 
discussion of the physical mechanism underlying this result is given in Sect. 4. 

12 



3.3 Comparison of (0 O ) an d (AA) 

In this paper we are not concerned so much with the numerical value of (0o) ; but 
with the fact that this condensate violates the SUSY indentity eq. (|5|). Since we are 
dealing with a matter of principle, it is worthwhile to give an alternative derivation 
of the existence of SUSY violations. 

Examining the various leading order supersymmetric results found in the liter- 
ature || H we observe that, typically, there are diagrams that contain a different 
number of insertions of the Higgs field. Each classical field provides a factor of v and, 
to match dimensions, the diagram must contain compensating powers of p. The rel- 
ative weight of different diagrams is therefore p-dependent, and the supersymmetric 
result depends crucially on performing the p-integration. 

In the calculation of (0o) the two diagrams Figs. 1 and 2 have a different p- 
dependence too. But now the p-integration gives rise to a SUSY violating result. 
Introducing to the dimensionless variable 

s = 4?rVp 2 , (21) 

we can rewrite eq. (|20|) as 

A 4 y 2 r 
(4>o) = m* — ^ 2 ^ h logy J dse- s s 2 (ds - c 2 ) . (22) 

The numerical values of the constants are 2c\ = c<i = 1. The reason why we have 
introduced them will be explained shortly. 

Let us now discuss the m^-dependence of (AA). The diagrams that contribute 
to the leading log are closely related to Figs. 1 and 2. We only have to replace the 
4>q tadpole (Fig. 3) by a new sub-graph in which the two A zero modes go directly 
to the external point (see Fig. 4). The other parts of the diagrams as well as the 
symmetry factors are unchanged. In the case of the diagram related to Fig. 1, this 
amounts to replacing the g 2 p 2 /(8h) factor in eq. (||) by 2. The coefficient of the 
leading log of (AA) takes the form 

A 4 lQy 2 f 

(AA) = m* — — — logy / dse s s(ds - c 2 ) . (23) 
v g° J 

Notice that the expression in parenthesis is the same as in eq. (p2] ). This expres- 
sion corresponds to the parts of Figs. 1 and 2 that are unchanged in the calculation 
of (AA). A comparison of eqs. (|22|) and (^) reveals that the integrand of eq. ( p2|) 
contains an extra power of s, i.e. an extra power of p 2 . Consequently, regardless of the 
numerical values of the constants C\ and C2, the leading SUSY violating logs cannot 
vanish simultaneously for (0 O ) and (AA). 
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As we have explained above, the power of p 2 in every graph is determined by the 
discrepancy between the explicit factors of v and the dimensionality of the diagram 
(taking into account the common dimensionful constant m^A 4 ). In the case of (AA), 
the integration over the pair of X s zero modes gives rise to no factors of v. In the 
case of the 0o tadpole, there is a factor of v 2 coming from the two insertions of the 
Higgs field indicated explicitly in Fig. 3. But this v 2 is cancelled by a 1/v 2 coming 
from the zero momentum Fourier transform of the 0o propagator: m~ 2 = (2h 2 v 2 )^ 1 . 
This non-analytic dependence on v 2 arises because the massive propagator is obtained 
from the massless one by an infinite sum over an arbitrary number of mass insertions. 
The overall power of v is eventually the same in both cases. Since the dimension of 
(0o) is smaller than the dimension of (AA) by two, there has to be an extra power of 
p 2 in the case of ((/> ). 

Substituting the numerical values of c\ and c 2 and performing the p-integration, 
we find that the coefficient of the leading log cancels between the two terms in eq. (|23|) . 
Thus, contrary to the case of (</>o), there is no logarithmic enhancement for (AA). 
A calculation of the contribution with no logarithmic enhancement to (AA) should 
involve the following ingredients. First, the cancellation of logarithmic IR divergences 
between the two diagrams discussed above leaves behind a finite non-zero remainder. 
Additional contributions arise from massless diagrams in which pairs of classical fields 
have been replaced by propagators. These diagrams are individually finite. 

All the diagrams mentioned above belong to the same order in the massless 
expansion. But this is not the end of the calculation. Comparable contributions 
arise from higher order terms in the Born series. Although diagrams with more 
insertions of the Higgs field are formally of higher order, they contain stronger IR 
divergences that compensate for the extra powers of the coupling constant y. We first 
observe that even terms in the Born series represent corrections to the ^-component 
of the "lepton" zero modes, whereas odd terms are corrections to the £- component. 
Denoting the (2n — l)-st correction by ^ n '(x) one has 

e {n) (x) ~ p mf" 1 ! 2 "- 4 , p < \x\ < m{ x . (24) 

In this equation we have shown only the leading power law behaviour and suppressed 
logarithmic factors. We see that the Born series represent an expansion in powers of 
m\X. But the physical cutoff of the x-integration occurs at x ~ m\ x . Thus, after 
the re-integration each extra power of x 2 leads to an extra power of 1/mf. A similar 
pattern is found for higher order corrections to the induced scalar fields (the sub- 
graph in the upper-left part of Fig. 2). The final result is that the massless expansion 
breaks down because it becomes an expansion in m\/m\ = 1. 
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The conclusion is that, in order to calculate the leading order value of (AA), 
one must sum contributions from all orders in the massless expansion. Within the 
massive Feynman rules, this is equivalent to a calculation of the non-logarithmic 
terms that arise from the exact form of the zero modes and the induced scalar fields 
in the "lepton" sector. We expect that eventually (AA) will be non-zero because the 
linearized massive field equations are manifestly not supersymmetric. 

4. Discussion 

The definition of the path integral measure requires one to specify a complete set 
of quantum modes. In the vacuum sector, the free (massive or massless) field equa- 
tions define the same set of plane waves for both bosons and fermions. Consequently, 
there are no SUSY violations in the vacuum sector. 

In the instanton sector, there are three relevant complete sets of modes. The 
first complete set consists of the massless modes pertaining to the self-dual classical 
instanton solution. This basis is universal for all spins, and it defines the formally 
supersymmetric perturbative expansion. The massive field equations, that contain 
the Higgs field, define two more complete sets: one for the bosons and one for the 
fermions. The massive bosonic and fermionic modes are different from each other, as 
well as from the massless modes. The path integral measure must be defined using the 
basis of massive modes, because the alternative basis of massless modes leads to IR 
divergences. 

In fact, the massless spectrum is not exactly supersymmetric. The massless 
perturbative expansion is well defined only in a finite box. Since there are no su- 
persymmetric boundary conditions in the instanton sector, the result is an 0(1/ R) 
discrepancy between the bosonic and fermionic spectra (R is the size of the box). The 
explicit 1/R breaking terms are multiplied by positive powers of R coming from the 
IR divergences, and so they can (and do) leave behind a finite effect in the infinite 
volume limit. 

The massive fluctuations spectrum, which is the appropriate one in the infinite 
volume limit, is manifestly not supersymmetric. The finite discrepancy between the 
spectra of massive fluctuations arises because of the non-constancy of the Higgs field 
in the instanton's core. The crucial effect is the existence of terms that involve 
derivatives of the classical Higgs field in the Schrodinger-like operators that define 
the bosonic and fermionic modes. (In the case of the fermions we refer to the square 
of the Dirac operator). Derivatives of the Higgs field appear as potentials in these 
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Schrodinger operators, and one can explicitly check that the potentials are different 
for bosons and fermions for any choice of the Higgs field. 

Let us now reexamine the SUSY violating result found in Sect. 3. The existence 
of different powers of p in the final expressions eqs. (|22f) and (|23"D means that the 
saddle point of the p-integration is different in each diagram. In other words, the 
effective value of p depends on the operator that we measure. This is not surprising. 
The integration over p is gaussian, and so p really represents a quantum mode. 

The logarithmic terms arise from the product of two wave functions that behave 
like 1/x 2 in the intermediate region. They depend on the effective value of p in two 
ways. The most important effect is that the normalization of the relevant wave func- 
tions is p-dependent. A secondary effect is that the short distance cutoff of the 1/x 2 
behaviour is determined by p. In the intermediate region, the equations of motion are 
locally supersymmetric to a first approximation. The SUSY violation arises because 
one cannot solve the equations with supersymmetric boundary conditions both a the 
origin and at infinity. Physically, the asymptotic region is supersymmetric because it 
is a free vacuum. The breaking of SUSY originates from the Higgs dependent part of 
instanton's core. In the calculation, it takes the form of non-supersymmetric values 
for the effective p. 

The Higgs model and the observables discussed in this paper were chosen for 
reasons of technical convenience. In essence, the anomalous breaking of SUSY was 
shown to be directly related to the IR divergences of massless perturbation theory. 
We thus conjecture that the SUSY anomaly is a generic property of asymptotically 
free supersymmetric gauge theories. 

Specifically, we proved the failure of the holomorphicity of physical observables 
constructed from the lowest components of gauge invariant chiral superfield. The 
conventional treatment of supersymmetric QCD depends heavily on the assumed 
holomorphicity, leading to conclusions such as the existence of a pathological run- 
away behaviour. We believe that these conclusions have to be reexamined in view of 
the failure of holomorphicity found in this paper. 
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Appendix 



I. Zero modes in the pure instanton background 

In the instanton sector, the classical gauge field is given by 



M 9 r 



and the classical Higgs field is 



< j, iA = iv^- ( p(r). (A.2) 
r 

We have suppressed the collective coordinate Xq. For mr <C 1 one has 

a(r) = r / \r 2 + p 2 ) , (A.3) 

(p(r) = r/{r 2 + p 2 )^ . (A.4) 

For mr ^> 1 these expressions are no longer valid. Both the gauge field and the Higgs 
field tend exponentially to their limiting values a(r) — > 1/r and <£>(r) — > 1. In this 
paper we only need the explicit form of the classical fields for mr <C 1. 

We next turn to the fermionic zero modes. The measure eq. (|7|) requires the 
normalization / <i 4 x|\I/| 2 = 1 for all zero modes. In the pure instanton background 
the normalized gaugino zero modes aref] 

{K)k ~ K(x 2 + p 2 ) 2 ' {A - b) 

where k = 1,2. The zero modes of the four charged doublets are identical except for 
the flavour quantum numbers. The Higgsino zero modes are 

(*PiAa)n = ~ 6 9 m6A ° 9 . 3 , (A.7) 



71" [X 2 + p 



where n = 1,2. The "lepton" zero modess are 



VAa± = W 2 , 2^ • (A ' 8) 
(x 2 + p 2 )2 



6 The euclidean partition function of supersymmetric theories is defined using a Majorana rep- 
resentation for the fermions. Consistency of this representation determines the global phase of the 
zero modes' wave functions. The Majorana representation is mandatory for the gaugino but it is 
convenient to treat the matter fermions analogously. See e.g. ref. for more details. 
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II. Exact zero modes 

When we turn on the Higgs field, the X ss and the r)± pairs remain exact zero 
mode, but with new field component and modified wave functions. The X ss pair now 
takes the form |§ 

(K)k = °lJ (A.9a) 
{$iAa)k = itiJAkg + ie f3a x fl x u a^ i ap k h (A. 9b) 

(lpOa)k = iSakP (A.9c) 

The hat denotes radial functions. The "lepton" zero modes are now given by 

VAa± = 5 aA u (A. 10a) 

£ iaT = ±e ia w (A. 10b) 

The quantum numbers of the different field components of these zero modes can be 
found in Table 2. The conserved angular momenta in the instanton sector are 

Kl = Sl + L1 + T a , 

K% = S^ + L a 2 + F a . (A. 11) 



channel 


5*1 5*2 T F L K x K 2 


A a 


\ 1 \ 




| | | 0,1 \ 


ipo 


\ | 


VA 


\ \ 


6 


\ \ 



Table 2: Quantum numbers of the field components of the exact fermionic zero 
modes. 

The radial functions are obtained by solving ordinary coupled differential equa- 
tions. The equations for the lepton zero modes are 

u' + 3au = —miipw (A. 12a) 

w = —rriiipu (A. 12b) 
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The prime denotes differentiation with respect to r. The functions a = a(r) and 
ip = if(r) are defined in eqs. ( |A.1| ) and ( |A.2| ). The equations for the X ss zero modes 
are 

f' + Aaf = -(ji/V2)ipg (A.13a) 

g' +{l/r-2a)g+{a-l/r)h 1 = -V^pupf (A.13b) 

p' = (m/Vtyytn (A.13c) 

^i + (3/r)Ai + 3(o-l/r)^ = V2rrvpp (A.13d) 

where 

h 1 = g + 2r 2 h. (A.14) 

The mass parameters /z, m and mi are defined in Sect. 2. Notice that the pairs (/, g) 
and (p, hi) diagonalize the mass operator at infinity. 

The Born approximation amounts to neglecting the off-diagonal term on the 
r.h.s. of eq. ( |A.12a| ). One solves eq. ( |A.12b| ) with u given by the pure instanton 
expression. (Compare eq. ( |A.10a| ) and eq. ( |A.8|) ). The result is eq. (0). Similarly, 



solving eqs. (|A.13b|) and ( A. 13d ) with / given by the pure instanton expression (see 
eq. ( |A.5| )) and with p = gives rise to eq. (|T2|). 

As discussed in Sect. 3, the calculation of the leading log arising from the "lepton" 
zero mode eq. fllTf ) requires a knowledge of the scale at which eq. (|14]) is no longer 
valid. The asymptotic behaviour of the wave function is a falling exponential 

u = uj = Co teVe^. (A.15) 



The transition from the power law behaviour given by eqs. ( A.8| ) and (|TJ) to the 



asymptotic behaviour eq. ( |A.15|) occurs at the distance scale m{ . This is a general 
property that can be derived as follows. Let us approximate the wave function by 
eqs. ( |A.8| ) and ( |HD for r < r , and by eq. ( |A.15| ) for r > r . We thus replace the 



smooth transition region by a sharp transition at r = tq. Requiring that the two 
components of the wave function be continuous across the transition point gives us 
two equations for the two unknowns r and c . The solution is r = Oim^ 1 ) and 
Co = 0(y). This information is sufficient to determine the leading log in eq. (|13|). 

The physical cutoff of the integral in eq. (|19D occurs at the same scale. This 
is derived by applying a similar analysis to the bosonic field equations. The same 
techniques can be used to obtain the qualitative structure of the transition region for 
the gaugino zero modes, but this information is not necessary for our calculations. 



19 



III. Approximate zero modes and the massive propagator 

The last issue that we discuss is the approximate zero modes originating from the 
Higgsino and X sc zero modes, and their relation to the propagators of the massive 
Feynman rules. Going to a partial waves basis, let us denote by Go the massless 
propagator with the quantum numbers of the Higgsino and X sc zero modes, i.e. 
K\ — and K 2 = \. The zero modes themselves will be denoted by the generic name 
^o- The massless radial propagator Go obeys the differential equation 

H G (r', r) = r" 3 5(r - r') - |tf > (*o| ■ (A.16) 



Notice the second term on the r.h.s. of eq. (|A.16|) that projects out the zero modes 



On the other hand, the massive propagator G with the same quantum numbers solves 
the equation 

(H + V) G(r', r) = r~ 3 5(r - r') . (A. 17) 

Now there is a simple delta function on the r.h.s. of the equation, because there are 
no exact zero modes with the given quantum numbers in the presence of the Higgs 
field. 

We now show that at short distances, the exact propagator G must have a "large" 
component which represents the approximate zero modes. This will justify the use 
of the massless Feynman rules for the Higgsino and X sc zero modes in Figs. 1 and 2. 
We first observe that the matrix element of the Higgs field 

E 1 = (* \V\* ) (A.18) 

between the zero modes is small (Explicitly E\ = e^gv/yl where the indices k and 
n count the superconformal and Higgsino zero modes respectively). If the zero modes 
have turned into finite energy bound states \I/ with E ~ Ex, the large component 
would simply be \^Sf) E~ x Now, the support of the zero modes is inside a potential 
well of radius p, and so for r m -1 it should make no difference whether the zero 
modes have turned into true bound states or into resonances. Let us denote the first 
Born correction by 

^ 1 = -G V^ . (A.19) 

One can check by direct substitution that for r m _1 , the massive propagator G is 
related to the massless propagator Go by 

G = G + |#o + *i> E^ 1 (#o + *i| + 0(g) . (A.20) 



In the instanton's core, the r.h.s. of eq. ( A.20|) contains an 0(1/ g) term which is an 
antisymmetric product of the original zero modes. 
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Within the massive Feynman rules, the first factor in eq. ([8]) is regarded as a part 
of fermionic determinant, 1.6. £IS Bj part of the measure. Keeping track of the relation 
between the diagrams drawn using the massive and the massless Feynman rules, we 
find that the application of the massless Feynman rules to the Higgsino and A zero 
modes in Figs. 1 and 2 is justified. 

Last we discuss the cancellation of logarithmic IR divergences related to the 
Higgsino and X sc zero modes. Such a logarithmic divergence can arise within the 
massless rules, if the two ip(x) wave functions which are to be integrated over at the 
ipip4 > Q vertex behave like 1/x 2 for p \x\. As an example, the integration over that 
vertex is logarithmically divergent in each of the three diagrams shown in Fig. 5. 
But if we first sum the three integrands (still using the massless Feynman rules) 
the integration over the vertex becomes IR convergent. (As in Sect. 3, the 

integration over that vertex is done after the integration over the external point 
eq. (p!l|). Notice also that the diagrams in Fig. 5 contain no logarithmic factors due 
to the vertex. The propagator lines, which replace the classical fields in Fig. 1, 
lead to a £(x) wave function which is damped at least by an extra 1/x 2 ). 

In Fig. 5(a), the logarithmic divergence arises from the first Born correction to 
the X sc zero mode, which behaves like ip(r) ~ 1/r 2 for r ^> p. Fig. 5(b) is a mixed 
case. In Fig. 5(c), the logarithmic divergence arises from the partial wave of the 
Go = (ifjift^ propagator with the same total angular momenta as the Higgsino and 
X sc zero modes. For r,r' ^> p, the propagator Go(r',r) has a piece that behaves 
like l/(r 2 r' 3 ). The ip end of the propagator behaves like 1/r 2 . (This is most easily 
seen by writing the propagator in a Dirac spinor basis). The reason for this unusual 
behaviour is the projection on the (Higgsino) zero modes on the r.h.s. of eq. ( A.16|) . 

Notice that the product of a Higgsino zero mode ip{r') times the first Born correc- 
tion ^(r) to a X sc zero mode behaves like l/(r 2 r' 3 ) too. An inspection of eq. ( |A.2(J| ) 
reveals that the massive radial propagator G(r', r) contains the two sources of a 
1 / (r 2 r /3 ) behaviour described above. 

We will now show that the two l/(r 2 r' 3 ) contributions on the r.h.s. of eq. ( A.2U| ) 
exactly cancel each other . The massive propagator G(r', r) is a solution of eq. ( |A.17|) . 
In the intermediate region rrT 1 ^> r, r' ^> p, one can neglect both the gauge field and 
the Higgs field. The dominant term in the equation is the free massless Dirac operator. 
Every piece of the radial propagator G(r', r) must therefore be the product of two 
homogeneous solutions of the free massless Dirac equation. For r > r', the possible 
terms are 1/r 3 and l/(r r') 3 . Corrections to these terms are damped by inverse powers 
of r and/or r'. In particular, a l/(r 2 r' 3 ) behaviour is inconsistent with eq. ( [A. 17f ) 
for r > r', and so it must cancel between the first and second terms on the r.h.s. of 
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eq. flPQP . 

Within the massive Feynman rules, the three diagrams of the massless Feynman 
rules (Fig. 5) are contained in a single diagram (Fig. 6). In this diagram one should 
use massive propagators. The thick lines now correspond to various components of 
exact zero modes. Sub-graphs describing the mixing of the Higgsino and X sc zero 
modes are absent, as they are regarded as a part of the the fermionic determinant. 
All other effects of the approximate zero modes are now contained in the massive 
fermion propagators. 

The cancellation of the l/(r 2 r' 3 ) terms on the r.h.s. of eq. ( |A.20| ) implies that 
the massive diagram Fig. 6 makes no logarithmic contribution to (0o)- We comment 
that, as a consequence of the absence of logarithmic terms, the sum of the three 
massless diagrams (Fig. 5) will be finite provided we amputate the external <p leg. 
Following the general pattern discussed in the introduction, that sum will be equal 
to the corresponding amputated massive diagram up to higher order corrections. 

As we have explained in Sect. 3, this cancellation occurs because subleading 
corrections to the wave function of an approximate zero mode cannot be disentangled 
from the continuous spectrum. This behaviour is natural if the zero mode has turned 
into a resonance. The new ip component of the \ sc zero mode is then responsible 
for its ultimate decay. In principle, there is also the possibility that the approximate 
zero mode has turned into a finite energy bound state. But since the energy of such 
a bound state is 0(m), it should behave like a resonance in the intermediate region 
m -i ^> r ^> p. The ip component of the X sc zero mode should then starts decaying 
outside of the potential barrier of radius p, but eventually it would remain trapped 
due to the mass term that changes the asymptotic value of the potential. 

The massless diagrams shown in Fig. 5 involve the same m^-dependent vertex 
as Fig. 1. There is a similar set of massless diagrams with the mQ-dependent vertex 
of Fig 2. Their sum is contained in the massive diagram shown in Fig. 7. (These 
massless diagrams can be retrieved from Fig. 7 in the same way that the massless 
diagrams of Fig. 5 are retrieved from Fig. 6). As before, the massive diagram Fig. 7 
makes no logarithmic contribution to (0o)- Finally, in Fig. 8 we show two massive 
diagrams that contain a closed loop. Their sum is UV finite, and neither diagram 
contains a logarithmic IR term. 
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Figure captions 



1. One of the two leading log contributions to (</> ). The diagrams in Fig. 1 up 
to Fig. 5 are drawn using the massless Feynman rules. The thick lines that 
emanate from the instanton (the shaded circle) represent fermionic zero modes. 
The dotted lines that end with a cross represent insertions of the classical Higgs 
field. The vertex marked with a thick cross is linear in TOq. 

2. The other leading log contribution to (4>o)- 

3. The 0o tadpole: the common lower part of Figs. 1 and 2. 

4. In the calculation of (AA) this sub-graph replaces the O tadpole in Figs. 1 and 2. 

5. Cancellation of logarithmic IR divergences related to the superconformal and 
Higgsino zero modes. In each diagram the integration over the ipip(f>Q vertex is 
logarithmically IR divergent. As shown in part III of the Appendix, the loga- 
rithmic IR divergences cancel among these diagrams before the p-integration. 

6. This diagram, drawn using the massive Feynman rules, corresponds to the sum 
of the three massless diagrams in Fig. 5. It makes no logarithmic contribution 

to (0 O ). 

7. Another massive diagram that makes no logarithmic contribution to (<f>o). This 
diagram contains the same m^-dependent vertex as Fig. 2. 

8. Two massive diagrams containing a closed loop. Their sum is UV finite, and 
neither diagram contains a logarithmic IR term. 
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